Thermodynamic coherence of the Variational Average- Atom in Quantum Plasmas 

(VAAQP) approach. 
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A new code called VAAQP (Variational Average-Atom in Quantum Plasmas) is reported. The 
model as well as main results of previous studies are briefly recalled. The code is based on a new 
fully variational model of dense plasmas at equilibrium with quantum treatment of all electrons. 
The code can calculate the Average Atom structure and the mean ionization from the variational 
equations respecting the virial theorem and without imposing the neutrality of the Wigner-Seitz 
sphere. The formula obtained for the electronic pressure is simple and does not require any numerical 
differentiation. A description of the principal features of the code is given. The thermodynamic 
consistency of the results obtained with VAAQP is shown by a comparison with another approach on 
the example of the aluminium 10 eV isotherm EOS curve. A first comparison to an iNFERNO-type 
model is also presented. 

PACS numbers: 52.25.Kn, 52.25.Jm, 52.27.Gr 



I. INTRODUCTION 



Models of plasma at local thermodynamic equilibrium 
taking into account bound electrons are being developed 
since long time ago. Such models first appeared in as- 
trophysics in order to calculate equation of state (EOS) 
and radiative properties of stellar objects. Calculations 
taking into account atomic electrons, allowing one to esti- 
mate the effective ion charge and photon absorption cross 
sections are essential for the modeling of stars. The abil- 
ity to account for bound electrons is important since it is 
well known that even relatively low Z elements are not to- 
tally ionized at temperatures typical for the central parts 
of the sun. More recently plasma atomic models taking 
into account bound electrons have also become impor- 
tant in the inertial confinement fusion (ICF) research, 
especially for the hydrodynamics modeling of ICF tar- 
gets. The role of the atomic physics models and codes 
in this modeling is to provide plasma EOS as well as 
transport coefficients like opacities and resistivity. With 
the development of the laboratory dense plasmas such 
models are also useful in plasma spectroscopy. Among 
the first models of atoms in plasmas was the well-known 
Thomas-Fermi (TF) spherical ion-in-cell model proposed 
by Feynman et al. 1} . This model based initially on ideas 
from the solid state physics extended to finite tempera- 
tures appeared to be thermodynamically consistent and 
successful in many practical applications. It is based on 
a variational principle, respects the virial theorem and 
provides a relatively simple expression for the electronic 
pressure. The plasma equilibrium can be constructed 
for each values of the three parameters : atomic number 
Z, ion density n.i and temperature T. The approach of 
Ref. [H is to model the whole bulk plasma through a self- 
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consistent calculation of a chosen central average atom 
treated within the quasi-classical TF statistical model. 
Non-central ions and electrons of the plasma may be 
viewed as a homogeneous jellium existing at the border 
and beyond the boundary of the central ion spherical 
volume which is electrically neutral. The radius of this 
sphere is the Wigner-Seitz radius Rws related to the ion 

1/3 

density by the formula Rws = (3/(47rni)) ' . The elec- 
tron density of the jellium ng related to the number of 
free electron per atom Z* = uq jni is not known a priori. 
It is found from the ion-in-cell TF model self-consistent 
field (SCF) calculation and equals to the electron density 
at the radius of the WS sphere. Similarly the electronic 
pressure in this model appears to be the pressure of the 
homogeneous electron gas at the atomic sphere bound- 
ary. The TF EOS is still being used in many ICF and 
astrophysical applications (see for example 0). Its quasi- 
classical character does not allows it to be used for the 
modeling of the radiative properties. 

For this reason with the appearance and progress in 
computers there have been several attempts towards 
quantum extensions of the average atom model. The 
most widely used and successful model has been pro- 
posed in 1972 by Rozsnyai (see Q). In Rozsnyai's ap- 
proach the free electrons have still been treated using 
the TF approximation. This is the principal reason why 
this model is not based on a variational approach and 
does not respect the virial theorem. The first entirely 
quantum model has been the Inferno model proposed 
in 1979 by Liberman (see 0,01) in which both bound and 
free electrons were treated quantum-mechanically. Liber- 
man has derived the self-consistent-field equations for the 
electrons from a variational principle however the aver- 
age ionization in his model is obtained from the neutrality 
of the WS sphere assumed in analogy to the TF ion-in- 
cell model. Both bound and free electron wave functions 
in Liberman's model are normalized in the whole space. 
Liberman's model does not respect the virial theorem 
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and does not have a simple expression for the thermody- 
namic pressure. Instead the pressure has to be calculated 
via a numerical derivation of the free energy per ion with 
respect to the ion density. One of the essential points in 
this model is the neglect of the Friedel oscillations [1, 0| 
beyond the WS radius. 

These oscillations have been taken into account in the 
model AJCI ("Atome dans le jeUium de charge imposee") 
proposed by F. Perrot In this model again, similarly 
as in Inferno, all wave functions are normalized in the 
whole space. In addition to this, following ideas from ear- 
lier studies in the condensed matter on impurities in jel- 
lium, F. Perrot has abandoned the approximation of the 
neutrality of the WS sphere. He has assumed only that 
the total SCF potential be local due to screening. This 
assumption is equivalent to the neutrality of all charges in 
the whole space. A serious problem which appears is how 
to determine the jellium electron density since the model 
in principle allows one to find a variational solution in 
the space of four parameters : Z, rii, T, and the average 
number of free electron per atom Z* . In such a way the 
AJCI model requires an additional equation allowing one 
to determine Z* . F. Perrot has proposed to use Z* from 
the TF ion-in-cell model i.e. from outside of the AJCI 
variational scheme. Due to this additional condition the 
AJCI model does not respect the virial theorem. 

First satisfactory answer to the question how to formu- 
late a fully variational approach to the quantum atom in 
jelHum has been proposed in [§, [l^, [Hj. The difference 
with respect to the previous approaches [1, H, [1] is the 
use of the cluster expansion for the two quantities : the 
plasma free energy per ion and the total number of elec- 
trons per ion. In both cases the zero order and the first 
order terms are retained from this expansion. As is shown 
in the cited references (see also Sec. 2), this systematic 
approach allows one to construct the variational system 
of equations respecting the virial theorem and to get a 
simple expression for the electronic pressure analogous to 
the one from As concerns the SCF equations for the 
bound and free electron density, our approach is similar 
to the one of the AJCI model. The difference with respect 
to this model is the presence of a supplementary equa- 
tion allowing one to determine Z* . This equation stems 
from the variational relation requiring that the total free 
energy per ion be stationary with respect to the average 
electron density rig. This is the only variational equa- 
tion in which both zero order and the first order terms of 
the cluster expansion are accounted for. This equation 
has been absent in the AJCI model since in this model 
as well as in the Inferno model the zero order terms 
in the cluster expansion have been disregarded. In the 
case of the Inferno model the requirement that the free 
energy be stationary with respect to no is replaced by 
the requirement that the WS sphere be neutral. In the 
case of the present model the WS sphere is in principle 
not neutral. Let us mention also that in the case of the 
quasi-classical Thomas-Fermi expressions for the free en- 
ergy and total number of electrons per ion, the approach 



from [T3| leads correctly to the TF ion-in-cell model of 
Feynman et al. 

This paper presents the main features of our VAAQP 
code which is based on the theory presented in Rcfs. 

and reports some preliminary results from this 
code. A first exploratory study (see [ill) ^^'^ already 
shown the feasibility of calculations with VAAQP in the 
Warm Dense Matter (WDM) regime. It appeared that 
significant effects of the variational treatment could be 
expected in the WDM regime. In this paper, we focus 
on the thermodynamic coherence of the results obtained 
with the code and on a first comparison with an In- 
FERNO-type model. In Sec. 2 the main features of the 
model will be summarized. The code will be described 
and some of the main results of the previous study will 
be briefly recalled in Sec. 3. In Sec. 4 we will present 
the dependence of the electronic pressure versus density 
at constant temperature equal to 10 eV in case of an 
aluminium plasma in order to prove the thermodynamic 
consistency of our model. Sec. 5 is devoted to a first 
comparison with an iNFERNO-type model. 



II. VARIATIONAL THEORY OF 
AVERAGE-ATOM IN JELLIUM 

In the cluster expansion the thermodynamic averages 
are expanded in the numbers of ions taken into account. 
As said before we first focus our attention of the free en- 
ergy. In the zero order ("no ions") all ions are smeared 
out into a neutralizing background and the free energy 
per unit volume /o('t.o; T) is that of a homogeneous elec- 
tron gas at the temperature T and of an electron density 
no which at this moment is considered as unknown : 



no 



n,;Z* 



(1) 



no is also related in the usual way to the unknown chem- 
ical potential hq. 

In the first order we consider one ion while all other 
ions participate in the surrounding jellium. The central 
ion will be interacting with the jellium. The cluster ex- 
pansion gives : 

=n,j d\ |/»"+^<=«»™ (X; n„ Z, T; r) - /o (no, T)} 

(2) 

The subtraction of the zero order free energy in the in- 
tegral assures the convergence of the term. In this way 
we define the following expression for the free energy per 
ion : 



Fq Fq (no, Ui, T) = 



F(X;n„Z,T) = fo + Ai^i 



(3) 



(4) 



on-i- jellium 



{X;ni,Z,T]r) - /q (no,T)| 

(5) 
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The first order contribution to the free energy per ion in 
our model consists of the three terms : 



AFf 



^pxc 



(6) 



AF° is the free energy of independent electrons 
(kinetic-entropic term) in a trial potential V{r) which 
we assume to be localized. 

AFf is the electrostatic interaction term taking into 
account all charge densities. The electron density is no 
more homogeneous and depends on the position r. For 
large values of the radius r this density tends to the value 
of jellium electron density of the zero order uq. The ion 
contribution to the charge density is due to the central 
nuclear point charge and to the charge density due to 
the non-central ions. We assume that this charge density 
is equal to no beyond a cavity of radius R and is equal 
to zero inside this cavity. In this way following the idea 
from Refs 0, |^ we introduce a cavity into which the 
non-central ion charge density does no enter. The only 
difference is that in our approach at this moment the 
cavity radius is still arbitrary. In such a way we get for 
the electrostatic interaction term, in atomic units: 

d^r\ {n{r) - noe{r ~ R))- 



AFf = - J 



Z 1 j- ^3^, j n{r') - no9{r' - R) 



AFf is the exchange correlation term: 



AFf 



d^r {fxc in{r)) " fxc ("-o)} 



(7) 



(8) 



The total free energy per ion shall be minimized with 
respect to the trial potential V{r) and the average free 
electron density uq taking into account the two following 
supplementary conditions: 



The total charge neutrality: 



Z 



(fr{{n{r)-no9{r-R))} = Q 



(9) 



The " ionization model" which we get from the two 
first terms of the cluster expansion applied to the 
total number of electrons per ion : 



Zrii — no + Hi / (fir {n{r) — no} 



(10) 



The substitution of Eq. [TO] into Eq. [§] leads immediately 
to the conclusion that the cavity radius R should be the 
WS radius : R = Rws- 

The explicit minimization with respect to the trial po- 
tential V{r) and the average free electron density no is 



thus applied to the thermodynamic potential 
n {no, Vir); n„ Z,T) = F {X; n„ Z, T) 



l[Z+ d^r{n{r)-noe{r~Rws)} 



(11) 



where 7 is a Lagrange multiplier. 

This minimization procedure leads, in the Local Den- 
sity Approximation (LDA) to the SCF equations: 

V{r) = V,i{r) - {n{r)) + V^, {no) (12) 



7 = ^0 + Vxcino) (13) 
where Vei{r) is the electrostatic part of the potential: 

,3 / f "(^') - na0{r' - R) ' 



and 



Va:c{n) 



dfxc{n) 
dn 



(14) 



(15) 



As said before, we get also a condition which allows 
one to determine no — n^Z*: 



d\{V,i{r)e{r~Rws)}^Q 



(16) 



Explicit analytical differentiation with respect to ion 
density of the thermodynamic potential at the equi- 
librium leads to the following formula for the electronic 
pressure: 



P = nj 



dni 



expl. 



= - fa + no ifJ-o + Vxc{no) + Vei{Rws)) 

(17) 



III. THE VAAQP CODE 

The VAAQP code solves the equations of the model 
described above. In the present version first it finds solu- 
tions of the SCF equations in the four-parameter space: 
{Z, T, ni, Z*) i.e. it proceeds as if an equilibrium could 
exist for any chosen value of the mean ionization Z*. 
Then it finds among these solutions the one that corre- 
sponds to the correct equilibrium i.e. such solution which 
fulfills Eq. [TH 

This code can calculate the electron density and 
the Average Atom structure using semiclassical (TF), 
quantum-non-relativistic (Schroedinger) or quantum- 
relativistic (Dirac) formalism. Finite difference methods 
are being used to compute electron density as function 
of the radius r up to an asymptotic radius at which the 
electron density is matched to its asymptotic form taken 
from the finite-temperature linear response theory [^, 0] . 
The iteration scheme previously used by Jena and Singwi 
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FIG. 1: Comparison between mean ionizations obtained from 
variational calculation (VAAQP), from non- variational NWS 
calculation (NWS), and from TFD calculation (TFD). 



(see [I4I) is used to obtain the self-consistency between 
the electron density and the potential V{r). 

At present, the code can use lyetomi and Ichimaru 
finite-temperature exchange-correlation free energy term 
[Tsi] as well as Dirac exchange energy term [13, [IE] but 
any other local-density exchange-correlation free energy 
term can also be used. In all the preliminary studies 
that we present in this paper, Dirac exchange term was 
chosen in order to provide direct comparisons with the 
Thomas-Fermi-Dirac (TFD) model. In the temperature 
regime of interest in these studies, the differences between 
results obtained using lyetomi exchange-correlation term 
and using Dirac exchange term appeared to be relatively 
small (the effect of finite temperatures seems to be small) . 

As stated before this code is primarily designed to cal- 
culate solutions of our variational model (which will be 
denoted by VAAQP). It has however other options. One 
of these allows one to perform the atom-in-jcUium cal- 
culation replacing the variational equation by the neu- 
trality of the WS sphere (the corresponding calculation 
will be denoted by NWS). Another option allows one to 
find iNFERNO-type solutions in which all calculations are 
limited to the neutral WS sphere. In this way, we can 
compare results from different approaches using exactly 
the same numerical procedures for the calculation of the 
electron density and other structure quantities. 

The VAAQP code was first presented in llj and it has 
been shown there using the example of the isochoric EOS 
curve corresponding to Al at 2.7 g.cm~^ matter density 
that variational treatment could lead to substantial dif- 
ferences of the calculated thermodynamic quantities with 
respect to results from non-variational models. These dif- 
ferences could be of the same order as the effect of the 
exchange-correlation or the effect of the quantum treat- 
ment of the electrons. They appeared to be especially 
important in the case of the WDM regime which is char- 
acterized in our model by weakly damped Friedel oscil- 
lations. As an example, we recall FIG. [T] which presents 
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FIG. 2: Comparison between electronic pressures obtained 
from Eq.[l7]and by numerical differentiation of the free energy 
along the 10 eV isotherm curve for variational calculations (a) 
and non-variational NWS calculations (b). 



the mean ionization as defined in Eq.[l] along the Al 2.7 
g.cm~^ isochoric curve for three models : VAAQP, NWS 
and TFD. Differences can be of the order of 30% and it 
can be seen that below 2.25 eV temperature, even the 
shell structure of the Average Atom is affected. 



IV. DIRECT VERIFICATION OF THE 
VARIATIONAL PRINCIPLE ALONG THE 
ALUMINIUM lOEV ISOTHERM CURVE 

The electronic pressure formula from Eq.[T7]is obtained 
by differentiation of the equilibrium free energy with re- 
spect to the ion density rij taking into account the con- 
ditions stated in Eqs. [TOl [51 As it is well known (see for 
instance [lo|). this is equivalent to the explicit analytical 
differentiation of the grand potential Vl with respect to 
the ion density. Each solution of a calculation with the 
variational option of VAAQP corresponds to the mini- 
mum free energy respecting Eqs. 1101 [9] i.e. the explicit 
minimum of O. Therefore, direct numerical differenti- 
ation of the free energy obtained with the variational 
option of VAAQP along an isotherm curve will also give 
the electronic pressure. In the case of a non-variational 
calculation such as with the NWS option of VAAQP, the 
grand grand potential f2 is not minimized. Calculations 
of the electronic pressure from Eq. [T7] and from numeri- 
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FIG. 3: Comparison between electronic pressures (a) and 
mean ionizations (b) obtained along the 10 eV isotherm curve 
for variational calculations (VAAQP) and non-variational 
NWS calculations (NWS). Pressures presented on (a) for 
NWS calculations stem from the numerical differentiation of 
the free energy. 
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FIG. 4: Comparison between electron density obtained for Al 
at 0.25 eV temperature and 2.7 g.cm~^ matter density from 
a variational calculation (VAAQP), a non-variational NWS 
calculations (NWS), and an Inferno- type calculation (IN- 
FERNO). Global plot of the electron density is presented, as 
well as enlarged views of the shell structure, Friedel oscilla- 
tions and vicinity of the WS radius. 



cal differentiation of the free energy should thus lead to 
values that are different. In case of a non-variational ap- 
proach none of these constitutes a justified way to obtain 
the electronic pressure as it is defined in thermodynamics. 
Nonetheless, numerical differentiation of the free energy 
is often used as a kind of a " mechanical" definition of the 
electronic pressure in some existing codes. 

In order to check the thermodynamic consistency of 
VAAQP variational option we can compare electronic 
pressures calculated from Eq. [IT] to the one obtained by 
numerical differentiation of the equilibrium free energy. 
Checking the agreement between these two quantities is 
not only a proof of the variational character of the model 
but can also be viewed as a very severe test of the perfor- 
mance of the code. We recall that such thermodynamic 
consistency has never been possible with other existing 
quantum Average Atom models. 

Comparisons between electronic pressures calculated 
from Eq. [17] and from numerical differentiation of the 
free energy along the Al 10 eV isotherm curve is dis- 
played on FIG. [U The case of VAAQP variational op- 
tion is addressed on FIG. [2^ whereas the case of the non- 
variational NWS option is presented on FIG. ^Bp. As 
expected, the agreement is found in the case of the vari- 
ational option whereas in the case of the non-variational 



option a disagreement can be observed. It is impor- 
tant to note that these differences tend to increase in 
the WDM regime. As it has been already noted in [ll[, 
rather important effects of variational treatment are to 
be expected in such regimes where the Friedel oscillations 
decay on a scale much larger than the WS radius. 

It is also interesting to compare the electronic pres- 
sures obtained from VAAQP variational option to those 
calculated by numerical differentiation of the free energy 
in the case of the non-variational NWS option. As it 
can be seen on FIG. [3^, these two electronic pressures 
obtained along the Al 10 eV isotherm curve are rather 
close. Nevertheless mean ionizations obtained from these 
two options can differ significantly, of the order of 10% 
(see FIG. |3}d). It means that average electron densities 
obtained from the two models can differ significantly de- 
spite the fact that both option lead to electronic pressures 
that are relatively close. It seems therefore that compar- 
isons between models should thus be made not only as 
concerns the EOS curves but also by comparison of other 
atomic physics properties involving the atomic structure 
such, for instance, as electrical conductivity. 
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V. FIRST COMPARISON WITH 
INFERNO-TYPE MODEL 

As a first step to comparison with iNFERNO-type mod- 
els, we present in FIG. 2] a comparison of the electron 
density obtained for the three models corresponding to 
the VAAQP variational option, the non- variational NWS 
option, and the Inferno option. 

Again, it is worth to note that results from the varia- 
tional option and from the non-variational NWS option 
differ strongly in the region dominated by Friedel oscil- 
lations. It is interesting to see that although Inferno 
model is based on the neutrality of the WS sphere, in 
the vicinity of the WS radius the iNFERNO-type electron 
density seems to be much closer to the variational solu- 
tion than to the electron density from the NWS option. 
The effect of the iNFERNO-type treatment seems to be 
important for the shell structure which appears to be 
different both from the ones of the VAAQP option and 
of the NWS option. 

VI. CONCLUSION 

In this paper, results of a new code based on a fully 
variational model of dense plasma at thermodynamical 
equilibrium are reported. This model takes into account 
the behavior of the self-consistent density and potential 
beyond the WS radius and allows one to obtain thermo- 
dynamically coherent value of the mean ionization. 



Thermodynamic consistency of the results obtained 
with the code is shown on the example of the aluminium 
10 eV isotherm curve. The results from this case show 
also that different models can lead to similar electronic 
pressures while providing significantly different average 
electron densities i.e. mean ionizations. 

The work on comparisons between thermodynamical 
quantities and atomic physics properties obtained from 
the variational VAAQP option and from the Inferno- 
type models is still in progress. As a first step, compar- 
ison of the electron density obtained with three models 
including an iNFERNO-type model is given. 
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